arXiv:1503.05247vl [astro-ph.HE] 17 Mar 2015 


To APPEAR IN ApJ 

Preprint typeset using L3T^ style emulateapj v. 5/2/11 


A SPECTROSCOPIC STUDY OF THE EXTREME BLACK WIDOW PSR J1311 -3430 

Roger W. Romani*, Alexei V. Filippenko^, and S. Bradley Cenko^’'* 

To appear in ApJ 

ABSTRACT 

We report on a series of spectroscopic observations of PSR J1311—3430, an extreme black-widow gamma- 
ray pulsar with a helium-star companion. In a previous study we estimated the neutron star mass as Mns = 
2.68 ± 0.14 Mq (statistical error), based on limited spectroscopy and a basic (direct heating) light curve model; 
however, much larger model-dependent systematics dominate the mass uncertainty. Our new spectroscopy 
reveals a range of complex source behavior. The variable He I companion wind emission lines can dominate 
broad-band photometry, especially in red hlters or near minimum brightness, and the wind flux should complete 
companion evaporation in a spin-down time. The heated companion face also undergoes dramatic flares, reach¬ 
ing ~ 40,000 K over ~ 20% of the star; this is likely powered by a magnetic held generated in the companion. 
The companion center-of-light radial velocity is now well measured with ATcoL = 615.4 ± 5.1kms“*^. We 
detect non-sinusoidal velocity components due to the heated face flux distribution. Using our spectra to excise 
flares and wind lines, we generate substantially improved light curves for companion continuum fitting. We 
show that the inferred inclination and neutron star mass, however, remain sensitive to the poorly constrained 
heating pattern. The neutron star’s mass, Mns, is likely less than the direct heating value and could range as 
low as 1.8 Mq for extreme equatorial heating concentration. While we cannot yet pin down Mns, our data 
imply that an intrabinary shock reprocesses the pulsar emission and heats the companion. Improved spectra 
and, especially, models that include such shock heating are needed for precise parameter measurement. 
Subject headings: gamma rays: stars — pulsars: general 


1. INTRODUCTION 

PSR J1311 -3430 ( hereafterJ1311)isa “black widow” mil¬ 
lisecond pulsar (MSP: IMetsch et aD2012^ with spin-down lu¬ 
minosity E = 5x 10'^'^ I 45 erg s“*^ (for a neutron star moment 
of inertia 10’^^ I 45 gcm^) in a Pb = 93.8 min orbit (iRomanil 
I 2 OI 2 I) . the shortest of any confirmed spin-powered pulsarQ 
J1311 is particularly interes ting, as photometric and spectral 
studies (iRomani et al.ll201^ hereafter R12) show that it has 
a strong, violently variable evaporative wind driven from a 
M 2 = 0.01 Mq He companion and that the neutron star is 
particularly massive. In R12 we found that with standard as¬ 
sumptions for heating of the companion by the pulsar, the 
light-curve modeling gave Mi = 2.68 ± 0.14 Mq (statis¬ 
tical errors), although with different assumptions for asym¬ 
metric heating one could obtain masses ^ 2.1-2.9Mq (sys¬ 
tematics range). This system promises to reveal much about 
close binary evolution ([Benvenuto. De Vito, & Horvathl2013l: 

Ivan Haften et al.ll2012l) . the efficacy of the companion evap¬ 

oration channel for producing isolated MSPs, and (possibly) 
fundamental constraints on the equation of state of matter at 
supernuclear densities. 

Measurement of the orbital motion, heating, and wind strip¬ 
ping of the companion is the key to such studies. To this end, 
we have made optical spectroscopic observations at several 
epochs. Here we summarize the measurements, describe the 
new phenomena discovered, introduce new models of the sys- 

* Department of Physics, Stanford University, Stanford, CA 94305- 
4060, USA; rwr@astro.stanford.edu 

^ Department of Astronomy, University of California, Berkeley, CA 
94720-3411, USA 

^ Astrophysics Science Division, NASA/Goddard Space Flight Center, 
MC 661, Greenbelt, MD 20771, USA 

^ Joint Space-Science Institute, University of Maryland, College Park, 
MD 20742, USA _ 

^ J1653—0158 IRomani et 3111201^ may be a similar system with even 
smaller = 74.8 min, but it has not yet had a pulse detection. 


tern, fit for the component masses, and comment on the re¬ 
sults. 

2. SPECTROSCOPIC CAMPAIGNS 

In R12, we discussed photometric evidence for substan¬ 
tial optical flaring activity of J1311. These flares appear 
to have a red spectral energy distribution (SED), but further 
study of this behavior requires spectroscopic characteriza¬ 
tion of the variable component. Also, multi-orbit photom¬ 
etry shows variable asymmetry about the optical maximum. 
This limits the reliability of model fitting (with simple sym¬ 
metric pulsar heating) in determining the system inclination i 
and the correction factor ATcor = ATcom/ATcoL, where ATcoL 
is the observed radial-velocity amplitude from spectral fea¬ 
tures on the heated face of the companion (the center of light; 
CoL) and ATcoM is the true radial-velocity amplitude of the 
companion center of mass (CoM). These factors are critical, 
since for a given observed radial-velocity amplitude the pul- 
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of the photospheric absorption lines can probe the surface- 
temperature distribution, and emission-line components can 
probe the evaporating wind. 

Thus, intensive spectroscopic campaigns offer the best hope 
for further understanding this important system. We report 
here on our efforts to date with the Keck-I/Low Resolu¬ 
tion Imaging Spec trometer (LRIS; O ke et al. 1995) and the 
Gemini/GMOS-S (iHook et al.ll2005h systems. For both, we 
sought long-slit multi-orbit coverage with high resolution and 
sensitivity. By placing the slit at a position angle (PA) of 
10°, we were able to simultaneously observe the pulsar and a 
r' = 18.05 mag G-type comparison star 15"fromJ1311 (Fig- 
ure[T]i. The latter allowed us to stably cross-calibrate the flux 
and wavelength, permitting coherent analysis of both datasets. 
With Keck-I we suffered from poor weather, with only half a 
night of the two scheduled nights delivering J1311 spectra, 
under marginal conditions. Nevertheless, we obtained three 
orbits of continuous coverage. The GMOS-S observations 
were queue scheduled and so all data were usable, with three 
visits covering slightly more than one orbit each. We describe 
these observations in detail, along with several short supple¬ 
mentary spectral sequences which give context to the variable 
behavior, especially near binary optical maximum. 

2.1. March 2013 LRIS Spectroscopy 

We observed J1311 with Keck-I/LRIS on 2013 March 
11 (MJD 56362.429-56362.626), using the l"-wide slit, a 
dichroic splitter at 5600 A, the 600/4000 grism in the blue 
camera, and the 400/8500 grating in the red camera. Table 
1 lists the wavelength coverage and resolution for all observa¬ 
tions. Our goal was to obtain high signal-to-noise ratio (S/N) 
spectra at all orbital phases. With a high southern declination, 
observations from Mauna Kea were perforce at moderate air- 
mass (secz = 1.7-2.5). The seeing was poor (mostly > 1.3") 
and the transparency variable. Each exposure was 300 s in du¬ 
ration; this was a compromise between improving the S/N but 
not having too much orbital smearing. 

Standard processing and optimal extraction were applied; 
fluxes were calibrated against the red-side and blue-side stan¬ 
dard stars BDh- 26°2606 and BDh- 28°4211, respectively. The 
comparison-star flux varied by a factor of 1.7 during this ob¬ 
servation sequence, largely owing to variable seeing-induced 
slit losses. After normalizing the flux, as monitored by the 
comparison star, we show (Figure |2l “trailed spectra” for the 
central portions of the blue and red sides. 

The 49 spectra cover just over three orbits ((j)B = —0.18 to 
2.85), with the midpoint of the barycentr ic arrival time of each 
exposure referenced to the ephemeris of iPletsch et al.l (l2012h . 
Here we set (ps = 0 at the first pulsar ascending node af- 
ter the start of obse rvations. For convenience, we will follow 
iBreton et al.l (120131) in referring to phase near optical maxi¬ 
mum (pulsar superior conjunction, (pB = 0.75) as “Day” and 
phases viewing the companion’s unheated face (pB = 0.25) 
as “Night.” “Dawn” and “Dusk” are at phases pB = 0.5 and 
pB = 0.0, respectively. 

Several features are evident in these trailed negative (darker 
means more flux) spectra. First, we see the continuum fad¬ 
ing and brightening from spectrum 1 (bottom) to 49 (top). 
The continuum modulation is stronger in the blue than in the 
red, as expected from the high temperature of the companion’s 
heated Day face. Superposed on this continuum are Doppler- 
shifted absorption features (R12). We also see broad and vari¬ 
able emission lines, strongest in the Dusk and Night phases. 
While little emission is visible during the first orbital mini¬ 


mum (spectra 5-10), the emission lines are very strong in the 
third Night (spectra 35^0), and indeed have such large equiv¬ 
alent width that they dominate in the red. Since the compan¬ 
ion is very H-poor (R12), both the photosphere and the wind 
trail driven from it are dominated by He I. 

2.2. 2013 Gemini Spectroscopy 

The Keck observations have good sensitivity and broad 
wavelength coverage, but lack the resolution to detect the line- 
shape distortions imposed by nonuniform companion heating. 
To search for such effects, we observed J1311 with Gemini, 
using the GMOS-S spectrograph in a novel mode to maxi¬ 
mize resolution in the blue part of the spectrum which dis¬ 
plays many photospheric He I and metal lines. These ob¬ 
servations used the R831 grating in second order through a 
0.75" slit, which covered 3950-5100 A at i? = A/5A 5000, 
^ 60 km resolution. This configuration required an order¬ 
blocking Alter. The best available Alter, g', unfortunately has 
a 4000 A cutoff which decreases sensitivity to Ca H; more¬ 
over, the protected silver coating on the primary mirror greatly 
limits its effective area below 4000 A. To minimize the or¬ 
bital velocity smearing, we restricted exposures to 175 s (i.e., 
maximum smearing of 60kms“^ at quadrature), obtaining 
38 consecutive exposures during each of three visits to the 
source. Each visit started near an optical maximum and cov¬ 
ered ApB = 1.27 in phase. The spectrum orbital phase is set 
by the exposure midpoint at the solar-system barycenter. 

The data were processed using the (Temini IRAE package. 
Eluxes were established using observations of the AO standard 
CD 32—9927. Comparison-lamp exposures were only taken 
once per visit, but, as for the LRIS spectroscopy, we moni¬ 
tored velocity stability and slit losses using the nearby com¬ 
parison star. This also facilitated comparison between observ¬ 
ing runs. The first visit, on 2013 May 6 (MJD 56418.15328- 
56418.23601), featured excellent ~ 0.5" seeing and high 
throughput, while the next two visits, 2013 June 3 and 5 (MJD 
56446.04081-56446.12281 and 56448.05900-56448.14101), 
had more typical ~ 0.8" image quality and lower S/N. 

In Figure |3] we show the average spectrum from maximum 
light (pB = 0.65-0.85, Day, all GMOS-S orbits) and near 
minimum light (pB = 0.15-0.35, Night, all GMOS-S or¬ 
bits), Doppler corrected to rest using the radial-velocity so¬ 
lution (see below). For comparison, the upper trace shows a 
He-dominate d model atmosphere with = 12,000 K and 
logp = 4.5 (iJefferv, Woolf, & Pollaccoll2001h . In the spec¬ 
trum from J1311 maximum, many narrow absorption features 
of He I and low-excitation metal lines are seen. In con¬ 
trast, the average spectrum at minimum shows broad emis¬ 
sion features, with He I and Mg I most prominent. Unlike the 
red LRIS spectra, these do not strongly dominate the 4000- 
5000 A continuum, so we cannot independently measure the 
wind emission-line profile and variation. In the trailed spec¬ 
tra this line emission is most obvious around pB ~ 0 (cen¬ 
tered on “Dusk” phases) during the first visit. It is weaker but 
present in the rest of the GMOS-S data. 

2.3. Other Spectroscopy 

Since the wind features sometimes show dramatic changes 
from one orbit to the next, we also describe several other 
observations of J1311, designed to check the system’s sta¬ 
tus. In R12, we discussed our initial six consecutive 300 s 
LRIS exposures covering pB = 0.56-0.93 on 2012 May 
17 (MJD 56064.306-56064.331). These data, obtained with 
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Fig. 2.— Trailed spectra from the 2013 March Keck/LRIS campaign. Fluxes have been normalized to remove variations in the simultaneously observed 
compaiison star. Upper panel: Blue side, showing Doppler-shifting absorption features from the companion photosphere. Lower panel: Red side, showing 


valuable He I emission features from the companion wind. The scale is inverted: 

the same configuration as on 2013 March 11, established the 
companion He dominance and provided an estimate Kohs = 
609.5 ± 7.5kms“^ for a sinusoid fitting the CoL radial veloc¬ 
ity. At this epoch J1311 appeared to be in quiescence, with no 
emission lines. 

We also observed J1311 on 2013 April 13 (MJD 
56390.424-56390.433) for 3 x 300 s near optical maximum 
((j)B = 0.75) using the DEep Imaging Multi-Object Spectro¬ 
graph (DEIMOS; Eaber et al. 2003) on Keck-II. DEIMOS 
provided coverage of 4450-9630 A, with the dichroic set at 
6990 A. Standard extraction, flux calibration, and cleaning 
were applied. Obs ervations were made at the parallactic angle 
(lEilippenkolll982h and did not include the comparison star. 

The companion was in a particularly bright state with a dra¬ 
matic flare during the second spectrum (centered at = 
0.835). During this flare the continuum flux increased by a 
factor of ^ 2. A He II absorption appears quite strong, par¬ 
ticularly in the second spectrum. Eigure|4]shows the dramatic 
continuum variability over ~ 300 s at this epoch. Note that 
the difference spectrum at bottom shows that the emission¬ 
line variations, if any, are quite small. We thus infer that this 
heating is applied below the photosphere, allowing relatively 
narrow absorption features to appear. The broader superim¬ 
posed emission lines arise from a different location. 

Einally, we returned to J1311 with Keck-I/LRIS on 2013 
May 10 (MJD 56422.347-56422.355) for three additional 
300 s exposures, this time covering phases (pB = 0.012- 
0.139. The source appears to have returned to quiescence with 


dai'ker indicates mpre flux. 

faint line emission. 

In Table 1 we summarize these various spectroscopic cam¬ 
paigns. To give some idea of the state of the wind line emis¬ 
sion, we list the flux of the He I A6678 line (we measure the 
He I A4471 line for GMOS-S spectra), after averaging over 
the full visit. Of course, as apparent from Eigure |2 the line 
strength is highly variable over even a single orbit, so these 
fluxes should be considered only a rough characterization of 
the wind activity. Eor example, although the line emission 
was not detected during the average of the second and third 
GMOS-S visits, we can detect He I emission during the Night 
{(pB ~ 0.25) portion of the orbit (Eigurej^l. 

Eigure |3 compares the XF\ spectra near (pB = 0.75 
from these five datasets. The GMOS observation and the 
LRIS/56064 spectrum (here and elsewhere, the label refers 
to MJD) have fluxes similar to that of LRIS/56362, but have 
been offset for readability. We also offset the LRIS/56422 
spectrum by —0.75; note that, unlike the other spectra in this 
figure, this is at quadrature {(pB = 0.01) and thus is intrin¬ 
sically fainter. Two DEIMOS/56390 spectra near maximum 
light are also shown: phase (pB = 0.77 before the main flare 
(only A > 5000 A plotted to avoid overlap) and the flare peak 
a.t(pB = 0.84. By (pB = 0.90, J1311 returns to the quiescence 
level. All quiescent spectra have similar continuum fluxes, but 
the He I lines vary greatly: absent on 56064, strong on 56362, 
and a factor of ~ 2 stronger on 56390. 

3. THE WIND KINEMATICS 
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TABLE 1 

PSR J1311-3430 Spectroscopic Campaigns 


Date 

(MJD) 

Tel./Instr. 

Range 

(A) 

Res 

(A) 

Exp. 

{N X s) 

Orb. Phase 
(4>b) 

(He I 6678) 

(10-18 gj.g g-l) 

56064^ 

Keck-I/LRIS 

3100-10,500 

4,7'’ 

6 X 300 

0.56-0.93 

< 3 

56362 

Keck-I/LRIS 

3100-10,500 

4,7'’ 

49 X 300 

-0.18-2.85 

200. 

56390 

Keck-II/DEIMOS 

4450-9630 

4.6 

3 X 300 

0.77-0.90 

390. 

56418 

Gemini/GMOS-S 

3950-5100 

0.9 

38 X 175 

0.58-1.85 

28.'’ 

56446 

Gemini/GMOS-S 

3950-5100 

0.9 

38 X 175 

-0.16-1.10 

< lO" 

56448 

Gemini/GMOS-S 

3950-5100 

0.9 

38 X 175 

-0.17-1.09 

< 8'’ 

56422 

Keck-I/LRIS 

3100-10,500 

4,7'’ 

3 X 300 

0.01-0.14 

9.0 


“ Previously described bv IRomani et alJ <201211 : all other obsei'vations presented here for the first time. 
Blue-side resolution. Red-side resolution. 

No red-side coverage; estimate from He I A4471. 


emission is intermittent along these vectors. 

Figure Q show a schematic model illustrating a possible 
kinematic origin for the wind line structure. Here we com¬ 
pute the cometary wind trail for matter stripped at low ve¬ 
locity (in the corotating frame) from the companion and then 
accelerated by pulsar radiation pressure; the only free pa¬ 
rameter in this model is opsr/^psr^ the ratio of wind ac¬ 
celeration to gravity. A similar model with the addition of 
a finite launch speed of 2/3 Vorb does a good job of defin¬ 
ing the radio ecl ipse envelope for the original black -widow 
PSR B19574-20 (Rasio, Shapiro. & Teukolsk^flO^ . Phase 
(fiB = 0.3 is shown in the upper panel, with the arrows and 
histogram indicating the possible radial velocities along this 
wind trail; a patch of He-containing gas will provide an emis¬ 
sion line within this envelope. The lower panel illustrates the 
range of velocities as a function of phase. This figure can be 
compared with Figure |6j we reproduce the two curves from 


a(A) _ 

Fig. 3.— Average spectra near maximum light {(f>B = 0.65-0.85, \ 
black curve) and minimum light (tj)B = 0.15-0.35, green curve) from the \ 
GMOS-S data. A comparison He-dominated model atmosphere (magenta "F 
curve) shows that many strong, narrow absorption features are present in § 
the maximum-light “Day” spectrum; some of the strongest lines are marked. "A 
Broad emission is apparent at minimum light. ^ 

The broad He I emission, which must arise in the compan- T 
ion wind, is clearly episodic and variable. Our best measure- ^ 
ment of this emission is in the March 2013 Keck/LRIS data 
(MJD 56362). Since the lines appear optically thin with rel¬ 
atively constant line ratios, we have formed a wind-velocity 
profile by taking the weighted, continuum-subtracted aver¬ 
age of the three strongest He I lines in the red; 5875.61 A, 
6678.15 A, and 7065.18 A. The result (Figure |6]l shows the 
flaring, varying emission-line profile. Scaled versions of this 
line profile provide good matches to the emission-line struc¬ 
ture in the blue, as well. 

In Figure |6] although the emission is patchy, there is ap¬ 
proximate orbital periodicity, with one component following 
the dotted sinusoid. Note that the emission other than this 
sinusoidal component seems to generally lie redward of the 
companion CoM velocity. We conclude that some preferred 
vectors of emission in the rotating system exist, but that the 



A(A) 

Fig. 4.— Sequential spectra from the blue side of the Keck-II/DEIMOS 
observations. He I emission from the wind is present in all spectra. The flux 
doubles during the second exposure (green curve), and the residual spectrum 
(cyan, at bottom) shows the absorption features after subtracting the contin¬ 
uum power law: He II absorption is strong. This flare component is only 
weakly present in the first spectrum (red curve), and no He II absorption is 
apparent during the third spectrum (blue curve) 600 s later. 
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that figure to guide the comparison. The wind trail plotted 
fades after A(j)B = 0.7; this boundary matches well in both 
amplitude and phase the quasi-sinusoid traced in Figure | 6 ]by 
the dotted curve. If the wind emission is relatively strong 
at lag A(j>B = 0.7, patchy at smaller A(j>B, but stronger at | 
(fiB ~ 0.75-0.25, this pattern provides a reasonable match to 5 o 
the data. Many orbits would need to be monitored to deter- | 
mine if this represents a dominant persistent emission pattern. 1 
This kinematic model is attractive, but the present data do “ 
not allow us to determine the physical origin of the emission¬ 
line velocities; of course, any simple sinusoid can be repro¬ 
duced by a rotating vector fixed in the binary frame, for ex¬ 
ample at the companion star. The dotted curve shown can 
be produced by a vector (i.e., a jet) from the companion with 
Vj ~ l.Svorb) in the star’s rest frame, directed in the orbital 
plane ^ 60° ahead of the line of centers. The physical nature 
of the pulsar flux driving the tail is also unclear. With E = 

5 X 10^^ I45 ergs“^, apsR/ffPSR = 3/4 requires a mean wind 
cross section of {ajm) Ri 1.5 x 10 ^cm^g“t. This is about 
right for a ionized pair plasma, but high for a baryonic wind. 

It is also interesting to consider what this wind tells us about 
pulsar evaporation. The momentum density of the pulsar wind 
is 40 /45 g cm“^ at the orbital radius. To ram-pressure 
balance this wind along the line of centers requires a compan¬ 
ion wind density pw = 1-0 x 10 “^^/ 45 t;//Q g cm“^ and a 
total mass flux rhw = ~ 6 x 10 ^^ / 45 t;/QQ gs“^, 

for a wind speed 600t;6ookms“^ comparable to the orbital 
speed and wind only from the heated face of the companion. 

Note that this gives an evaporation timescale r Ri TO 2 /tow ~ 
10 m_ 2 U 6 oo/-f 45 Gyr, for m 2 = lO“^m_ 2 M 0 . Although 
m 2 was larger when accretion stopped and evaporation be- 



Fig. 5.— Spectra near optical maximum {(^b = 0.65-0.85) from our 
five observing campaigns (grouping the three GMOS-S visits). The GMOS- 
S (cyan) and LRIS/56064 (blue) spectra have very similar flux to that of 
our long LRIS/56362 campaign (black) and are offset for readability. The 
LRIS/56422 (magenta) is also offset by —0.75, but is at “Dusk” phase and 
so is intrinsically fainter. The other spectra have flux as plotted. Note the 
dramatic continuum increase during the DEIMOS/56390 flare (green), dur¬ 
ing which He II absorption is prominent. Also note the strong and vaiiable 
He I emission in the red. 



Orbital Phase 


Fig. 6.— Velocity profile of the He I emission from the companion wind 
on MJD 56362. For reference, the red curve shows the radial velocity of the 
simple sinusoidal fit to the absorption features, marking the CoL velocity of 
the heated surface of the companion star. The emission lines do not track this 
motion, but a likely sinusoidal component is shown by the dotted line. Over 
three periods are shown^ g^ith dai'k colors ](i^ai'king bright lin^emission. 



0 0.5 1 1.5 2 

Orbital phase 0g 

Fig. 7.— A simple ballistic model for the cometaiy trail emission. Mass is 
lost from the companion and accelerated radially at 3/4 gpsR, via the pulsar 
spin-down flux. One phase, with the trail (out to a lag A0)s = 0.7), along 
with the line profile (histogram) for uniform illumination, is shown at top. 
The red arrows show the companion velocity in space and at the red edge of 
the line profile, for (ps = 0.3. Two periods of the trail velocity range along 
with the companion velocity and the sinusoidal component curves (curves 
exactly-as in Fisure[^ are shown in the lower panel. 

gan, E was Higher in the past so the companion may be com¬ 
pletely ablated in less than a spin-down time, especially for a 
slow wind or > 1 . 

4. THE PHOTOSPHERE FLARE VARIABILITY 

Figures (Hand [5] show that in addition to the emission-line 
variability, there can be epochs of true con tinuum increase . 
The dramatic i' flare shown in Figure 1 of [Roman] (l2Q12h . 
which erupts to 6 times the normal flux at 0 . 8 , is likely a 
similar event. Figure[4|gives an important clue to the nature of 
these variations. The difference spectrum shows that the flare 
does not affect the emission-line fluxes, but produces strong 
absorption lines. The He II absorption dramatically increases 
(including a very strong, unusual detection of He II A6560, 
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generally confused with the absent Ha line). We infer that 
these flares must be caused by impulsive deep heating well 
below the photosphere. However, since the photometric flares 
and spectral flares have decay times of ~ 300 s, if we assume 
a sound speed ^ 8 T 4 in the ~ 30,000K He gas seen 

during these bursts, we infer a characteristic depth of ~ 2 x 
10®cm(-0.03i?*). 

By measuring the He II/He I line-intensity ratios in the 
flare spectrum, compared with th e ratios in model spectra of 
iJefferv, Woolf, & Pollaccol (l200lb . we find an effective flare 
temperature of Tpi = 39,000 ± lOOOK. The continuum flux 
increase from quiescence is a factor of 2.2 at 5500 A assuming 
Planckian continua in this range. This implies a flare-emitting 
area 0.27 times that of the quiescent star. Thus, we infer that 
a large fraction of the star brightened by a factor of 8.3 during 
the eruption. At the estimated (from the dispersion measure) 
distance of 1.4kpc, the equivalent radius of the flare area is 
0.05 R 0 , which is 0.5 7?*, in good agreement with the spec¬ 
tral estimate. 

The average thermal luminosity of this burst is 5 x 
10^^ergs“^. This is ^ 80 times the pulsar flux hitting the 
companion, so we infer that these outbursts represent stored 
energy released deep below the photosphere. The natural can¬ 
didate energy source is a magnetic held, with characteristic 
magnitude (B) Ki 10 kG required to power the flare. The 
origin of such a held is unclear. Typically stars with T^ff = 
12,000 K are radiative throughout. However, here we expect 
that high-energy particles from the pulsar are deeply heating 
the photosphere, with (for example) TeV depositing their 
energy at S « 300gcm“^. While this does not in itself in¬ 
duce convection, we expect latitudinal flow because of the 
large temperature gradients from the differential heating. If 
convective, given the short effective spin period for the tidally 
locked companion, we have the ingredients needed for a ro¬ 
bust dynamo, large magnetic fields, and energetic flare events. 

5. IMPROVED QUIESCENCE LIGHT CURVES EROM SPECTRA 

The large equivalent width of the variable He I lines can 
certainly affect the broad-band fluxes, especially at binary 
minimum and for red colors. This is likely the origin of the 
red “flaring events” noted by R12 in simultaneous multi-band 
GROND light curves. We can use our spectra to partly cor¬ 
rect for this effect, synthesizing light curves from the contin¬ 
uum. To do this, we integrate our normalized spectra over the 
SDSS g'r'i'z' Alter passbands, after excising ±1500km s“^ 
wide bands around the strong wind lines (mostly He I, but also 
Ca and Mg) seen at binary minimum (Figure | 6 ]l. To normal¬ 
ize these Alter bandpasses to standard SDSS magnitudes, we 
integrated the comparison-star spectra over the same bands. 
Carefully calibrated GROND observations give magnitudes 
for this star of g' = 18.58, r' = 18.05, i' = 17.79, and 
z' = 17.70 mag (A. Rau, private communication). This de¬ 
termines an offset to the measured magnitudes for each spec¬ 
trum. Comparing with synthesized magnitudes for the full 
SDSS passbands, we do indeed And decreased source fluxes, 
by as much as a magnitude when the wind lines are strong. 

However, even this line excision does not completely re¬ 
move the light-curve fluctuations. To best estimate a “quies¬ 
cence” light curve, we use our spectra from MJD 56362, and 
compare the synthesized magnitude estimates in At); = 0.05 
phase bins. We adopt the largest (faintest) magnitude in each 
band in each bin as the quiescence value and plot the resulting 
light curves in Figure 0 Statistical errors are obtained from 
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Fig. 8.— Synthesized band magnitudes from the Keck MJD 56362 spectra. 
These g'r'i'z' magnitude estimates are formed by integrating over the SDSS 
passbands after excising strong (mostly He I) emission lines. Magnitudes are 
normalized by matching the synthesized comparison-star magnitudes to the 
direct imaging values. Two periods are shown with estimated errors during 
the second cycle (see text). 

the fluctuations in the spectroscopic fluxes integrated to derive 
the band magnitudes. We further estimate the remnant light 
curve fluctuations associated with faint wind lines or contin¬ 
uum flares, by taking the difference to the second faintest 
magnitude observed for a given bin. These error estimates 
are added in quadrature to the statistical errors and plotted on 
the bin points in Figure | 8 ] 

Several features of the derived light curves deserve com¬ 
ment. First, our synthesized magnitudes, normalized to the 
comparison star, are roughly 0.6 mag fainter than those re¬ 
ported from direct GROND imaging in R12. It is possible 
that this apparent 40% flux decrease is caused by differen¬ 
tial slit losses between the pulsar and comparison star. How¬ 
ever, the GROND magnitudes at MJD 56118 may have been 
well above quiescence, since these data show larger fluctu¬ 
ations (even at Ri 0.75) than seen in these synthesized 
magnitudes. Of course, even these new light curves will not 
represent complete quiescence — fluctuations are still visi¬ 
ble across the maximum — but they provide a substantially 
improved target for quiescence modeling. Overall the light 
curves follow the form of the direct imaging SOAR/GROND 
curves reported by R12, although we appear to have better 
captured the red colors at binary minimum. It should be noted 
that as in R 12 , the light curves remain slightly asymmetric, 
brighter at “Dusk” {(pB ~ 1.0) by ^ 0.1 mag compared to 
Dawn {(pB ~ 0.5). 


6. THE COMPANION RADIAL VELOCITY 

The dramatic heating offsets the CoL in a poorly un¬ 
derstood, temperature-dependent way. Also the variable 
wind emission in precisely the He I lines dominating the 
companion photosphere can shift the absorption-line cen¬ 
troids. These effects make measurement of the companion 
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radial velocity challenging. We focus here on the March 
2013 Keck/LRJS and the May-June 2013 Gemini/GMOS-S 
datasets. For the Keck data we start by removing a scaled 
version of the emission-line velocity profiles (Figure |6]l from 
each detected line feature in the blue channel. We then use 
a He-dominated model atmospher e template from the atlas of 
[Jeffery. Woolf. & Pollaccol (|20^ (T.ff = 12,000 ± 1000 K, 
log g = 4.6 ± 0.2) as a template for velocity cross-correlation 
to obtain the best radial velocity of each individual spectrum. 
Cros s-correlation is perfor med using the IRAF RVSAO pack¬ 
age (iKurtz & MinklTTOQ^ . The quality of the fit is moni¬ 
tored by the correlation statistic R (the ratio of the correlation 
peak height to random peak height), which varies from ~ 12 
on the day side, where velocity uncertainties are as small as 
ay = 15kms“^, to an insignificant i? < 1 on the Night side. 
Measurement of the comparison star with a G-star template 
showed a slow drift in radial velocity from H-30kms“^ to 
—30kms”^, with typical uncertainties of 6kms“^. We ac¬ 
counted for this residual calibration drift and also corrected to 
heliocentric velocities. 

Similar measurements were made of the GMOS-S spec¬ 
tra. Here we reach i? « 7 during the Day phases. With the 
high resolution, this delivered velocity uncertainties as small 
as ~ 5kms“^. Night-time cross correlations were of¬ 
ten insignificant, but occasionally delivered apparently sig¬ 
nificant R ss 1-2; some of these measurements were close 
to the orbital curve while others had large departures. The 
comparison-star radial velocity was monitored to typically 
ay ~ 4-6kms“^, although uncertainties grew to as large as 
~ 10kms“^ when conditions were poor. No strong velocity 
trend was evident in these data. 

The corrected cross-correlation velocities for these two data 



Fig. 9.— J1311 radial-velocity measurements from LRIS and GMOS-S 
data. The first cycle shows all > 1 coiTelation points, with velocity error 
Hags. During the second cycle the data used in the simple sinusoidal fit (solid 
line, dashed line marks the systemic velocity F) are in green and cyan; red 
points are excluded. The bottom panel show the residuals to this best-fit 
sinusoid. The cuiwe plots residuals to the ELC velocity curve (“ELC-Sp” 
model) after subtracting the simple sinusoid, demonstrating that the model’s 
non-sinusoidal radial-velocity terms are present in the data (see text). 


sets are show in Figure |9] Two periods are illus trated, phased 
to the orbital solution of iPletsch et alJ (l2012h . During the 
first, all points with i? > 1 are plotted, with their (t„ error 
flags. The radial-velocity curve during the Day phase is very 
well measured. During the Night phase a number of GMOS- 
S points appear to follow the expected curve, but many are 
displaced by 4-30-60km s“^; the few plotted LRIS points 
in this phase range seem largely random. Recalling that the 
GMOS-S data measure only the blue portion of the spectrum, 
it may be that the radial velocity is tracking flux from the 
bright lune of the heated side, visible at the companion limb 
for i < 90°. 

In addition to the low-statistics Night-time points, it is clear 
that there are large variations from the radial-velocity curve 
at the Dawn and, especially. Dusk phases. At times the ve¬ 
locities follow the simple sinusoid, at others the velocities lie 
up to 100kms“^ “inward” of the expected curve, to the 
blue (—v) near (j)B = 0.5 and to the red (-t-r^) near (j)B = 0.0. 
These offsets are highly significant and do correlate among 
spectra. However, we are not able to associate these with, say, 
periods of strong wind emission lines. One interpretation is 
that these offsets represent radial outflow near the compan¬ 
ion terminator; an emission component to the red (blue) at 
Dawn (Dusk) would impose a variable blueshift (redshift) on 
the absorption features measured against the cross-correlation 
template. More high-quality data are needed to pin down the 
physical origin of this effect. 

One consequence of the terminator shifts is a difficulty 
in fitting a simple radial-velocity curve. If we fit only 
to correlation measurements with R > 3, and restrict to 
(pB = 0.5-1.0, we obtain an observed radial-velocity am¬ 
plitude ATcoL = 611.4 ± 5.4kms“^ and systemic velocity 
F = 39.4 ± 3.7kms“^ (statistical errors only). The fit de¬ 
partures remain substantial, with a per degree of freedom 
(DoF) of 9.8 for 74 DoFs. If we restrict further to the Day side 
(0.55 < (/)B < 0.95), we obtain ATcoL = 615.4± 5.1kms“^, 
with x^/DoF = 4.0. In contrast, fitting all points with R > 3 
includes many discrepant Dusk points and gives ATcoL = 
590.7 ± 7.2kms“^, with x^/DoF = 15.6. It is evident that 
non-sinusoidal terms are present. The bottom of Figure |9] 
shows the residuals to the simple sinusoidal fit. Despite the 
scatter, some systematic trend is visible, especially during the 
Day side. The curve plotted, for a best-fit ELC model with 
an LI cool region (see below), shows a similar residual to the 
fit sinusoid, and the fit to this curve decreases x^ by 53.7, i.e. 
X^/DoF decreases from 9.8 to 9.1. 

7. MODELS OF THE HEATED PHOTOSPHERE 

There are a number of codes used to model light curves and 
radial velocities in interacting binaries. We And that none is 
fully adequate to cover the extremes of temperature, surface 
composition, and pulsar heating found in J131L Neverthe¬ 
less, such modeling fails in instructive ways and provides a 
guide toward future, more detailed solutions. We report here 
on modeling with the ELC code (Grosz & Hauschildt 2000) 
and the ICARLfS code (Breton et al. 2011), both kindly made 
available to us by their authors. Eor both we can run in “MSP” 
mode, with the period Pb = 0.065115d and pulsar orbit 
xi = aisini = 0.0105811t-s held fixed. Both codes model 
the effect of pulsar heating of the companion as illumination 
by a point (X-ray) source. 

Eor the ICARUS code we ran with color tables generated 
from the BTSettl atmospheres (Allard et al. 2010). In all 
cases the best-fit models preferred large Roche-lobe All fac- 
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Fig. 10.— Residuals to the model fits. Top: photometric {g'r'i') residuals 
with (left to right) the ICARUS fit, the ELC baseline fit, and the EEC fit with 
a 0.6 X Li spot. Bottom left: the radial-velocity points used in the ELC 
fitting. Middle and right panels show the residuals to the ELC baseline and 
Li spot fits, respectively. 

tors; we set the effective fill factor to 0.99. In fitting our 
quiescent g'r'i'z' light curves, the main parameters are the 
orbital inclination i, the Night side (unheated) star temper¬ 
ature T^v, and the heated side (Day) temperature Tq. We 
report this last parameter as an effective heating luminosity 
Lh = 47ra^cr(T^ — T^), where the orbital separation is 
a Ki 0.95 Rq and the heating luminosity can be compared 
with the pulsar spin-down power E = 5 x 10^"* I 45 ergs“^. 
ICARUS also fits the absolute fluxes and so returns a dis¬ 
tance estimate and an effective extinction. If left free, the 
extinction hts to an insignificantly small value; to be con¬ 
servative, we set the value to the low net Galactic extinction 
Av = 0.173 mag in this direction (Schlafly & Finkbeiner 
2011). The fit values are given in Table 2, with uncertainties 
as reported by the ICARUS fitter. The distance is surprisingly 
large at 2.6 kpc, and the heating flux slightly exceeds the stan¬ 
dard pulsar spin-down power. Most importantly, the ht incli¬ 
nation is so small that with the observed radial-velocity ampli¬ 
tude KcoL = 615kms“^, the minimum pulsar mass would 
be 2.77 M©. Since the true kinematic amplitude is larger (typ¬ 
ically with ATcor > 1.03, TablelU, then ATcoM > 633kms“^ 
and the inferred mass would be unphysically large at > 3 Mq. 
Also, the ht residuals (Figure fTOl i demonstrate that the max¬ 
imum is substantially wider than the model light curve, and 
that the colors are redder at peak than predicted by the model. 
Clearly, the present ICARUS model is inadequate for a de¬ 
tailed ht. 

The ELC code was used to ht the light curves and spectra 
simultaneously. For these hts we used an atmosphere model 
table computed for SDSS hlters by J. Orosz, extracted from 
the “NextGen” atmosphere grid wi t h an e xtension to Tgff > 
10,000 K using ICastelli & Kurucd (l2004ft models. As be¬ 
fore, the hts preferred large Roche-lobe hll factors and this 
value was hxed at 0.99. We again ht for the inclination. 


Tjv, and heating, but additionally constrained the mass ratio 
q = Ml/M 2 for direct estimates of the component masses. 
The results, which include radial-velocity constraints, prefer 
larger inclinations and lower masses than the ICARUS solu¬ 
tion. The Tjv temperature is slightly higher, but most impor¬ 
tantly the ht Lh is much higher and is indeed substantially 
larger than the total spin-down power. Again, the observed 
maximum is appreciably wider than the heating model sug¬ 
gests, and the mid-Day {(/ = 0.25) temperature is cooler than 
predicted. The ht is not statistically satisfactory and implies a 
large pulsar mass mi = 2.6!to 2 - Here, and for all ELC hts, 
we quote “Icr” errors as the full multi-dimensional projection 
over all ht parameters of the region about the ht minimum 
with a ht statistic increase of Xmm/HoF = 1; in this sense 
our uncertainties are conservative, with marginalization over 
all other parameters and error inhation for imperfections of 
the model. 

The ELC code has the option of applying artihcially heated 
or cooled star spots to the secondary. As in R12, we can sub¬ 
stantially improve the ht by cooling the inner Lagrange point. 
We apply a spot which decreases the local unperturbed surface 
temperature with a Gaussian prohle of radius 45°. The best 
hts invoke a ~ 40% temperature decrease and do give a sub¬ 
stantially better model for the light-curve maximum. The re¬ 
sult is that the companion is brightest in a ring pulsar-ward of 
the terminator, with the CoL much closer to the CoM than for 
direct heating (small ATcor)- In addition, the heating pattern 
broadens the light-curve maximum. This drives the inclina¬ 
tion i to larger values, to maintain the depth of the light-curve 
minimum. Both of these effects serve to decrease the inferred 
pulsar mass. 

While there are physical effects (see below) which can 
accommodate a larger effective heating power, Lh = 
1035.6 gj-g g-i jg larger than the pulsar can supply. Also, 
the large Tjv gives a poor match to the relatively red colors and 
large magnitude at c/b = 0.25. The drive to these large values 
is a consequence of the spot approximation for the surface- 
temperature distribution. If we hx Lh = 10^®ergs“^ and 
Tm = 4500 K, the largest values consistent with the observa¬ 
tions, we hnd that the inferred inclination (and thus neutron- 
star mass) remain a very strong function of the spot temper¬ 
ature decrement (Eigure [IB, which is rather poorly deter¬ 
mined. The best fit, with a 0.6 x spot (40% spot decrement), 
implies Mi = 1.9 Mq. Weaker spots lead to large inferred 
masses, stronger spots requir e edge-on orbit s in conflict with 
the lack of an X-ray eclipse (lRomanil|2012h . Thus, our con¬ 
clusions depend sensitively on details of the heating model. 
Notice in EigurefTOlthat the radial-velocity residuals vary only 
weakly between quite different heating patterns. The main ef¬ 
fect is an overall increase in amplitude, which is reflected in 
the changed ATcor; the temperature (color) is more sensitive to 
the heating model. 

8. CONCLUSIONS: INDIRECT HEATING, NEUTRON-STAR MASS 

Our goal of measuring the orbital inclination and the ATcor 
factor with light-curve fitting has been frustrated by the poor 
fit to the light-curve maximum. A strong clue to the nature of 
the difficulty is that simple direct-heating fits require a heating 
power larger than can be supplied by the 5 X 10^"^/45 ergs ^ 
spin-down power. Eor the ICARUS model. Lx exceeds the 
nominal spin-down power by a factor of 1.3. Eor the basic 
ELC model, the excess is a factor of 6, and for the ELC spot 
model, the heating power exceeds that supplied by the pulsar 
by a factor of 8. This is especially worrisome as the typical 
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TABLE 2 

J1311-3430 Light Curve/spectral Fits 


Parameter 

Units 

ICARUS 

ELC 

ELC-Sp 

ELC-Sp/ig 

Inch 

Deg. 

57.1 ± 1.5 

64.0 ±3.5 

8llf 

77tl 

Tn 

K 

4240 ± 150 

4800 ± 400 

5000 ± 400 

4500'= 

Log(Lif) 

erg s~^ 

34.80 ± 0.03 

QC C + O.IS 

00.0_Q 3 

35.6 ±0.2 

35.0'= 

Q 


- 

184.5 ± 1 

179.51J.S 

180 ±0.5 

^cor 


- 

1.062 

1.033 

1.036 


Mq 

- 


1 Q9+0-15 

±.o^_0 06 

-L-o»_o 06 


10-2 M0 

- 

1.42±0.15 

1 

-‘-•'^-*--0.02 

^•^^-0.03 

X^/DoF 

148 

- 

3.74 

2.60 

2.66 


“ Value infen'ed from fit - projected uncertainties, errors inflated by x^/DoF. 
Bounded by lack of an eclipse. 

Fixed in fit, with spot depth 0.6X. 


heating efficiency estimated by iBreton et al.l (1201 3h is 77 = 
0.15 of the pulsar spin-down. 

This discrepancy is so large that a number of factors are 
likely relevant. First, the effective moment of inertia and re¬ 
sulting spin-down power may be larger than assumed, espe¬ 
cially if Ml is > 2 MfTi where I/ih ~ 2 ^, d epending on 
the equation of state (iLattimer & Schut3l2005h . Also, the 
present modeling assumes normal composition in the atmo¬ 
sphere; He-dominated spectra can be bluer, allowing a lower 
fit Teff. Most importantly, the assumption of direct, radia¬ 
tive heating in the ICARUS and ELC models is likely inade¬ 
quate for J1311. Instead, we expect that the pulsar wind ter¬ 
minates at an intrabinary shock wrapped around the Li point, 
reprocessing the pulsar wind into high-energy particles which 
deliver heat below the companion photosphere. This repro¬ 
cessing can both increase the effective heating flux over di¬ 
rect radiation, as the intrabinary shock subtends a larger angle 
at the pulsar, and preferentially heats the companion termi¬ 
nator. Further, Coriolis forces break the shock symmetry in 
the rotating frame and may thus provide a mechanism for ex¬ 
plaining the asymmetric light-curve maximum. The presence 
of strong companion fields at the Li point (indicated by the 
large flares) may also preferentially direct intrabinary shock 
particles away from Li. Such effects are dramatic for J1311, 
but may be present at lower levels in other evaporating pulsar 
binaries, where hi gh-precision photometry sho ws asymmetric 
light curves (e.g.. lSchroeder & Halpernll2014l) . This may bias 
mass estimates inferred from fitting light curves with simple 
direct-heating models. 

Our spectroscopic campaign has revealed a number of new 
aspects of this remarkable system. First, J1311 undergoes 
violent flaring in the optical, which suggest strong mag¬ 
netic activity on the heated face of the companion. Second, 
strong (but intermittent) line emission is generated in the wind 
driven off the companion. This outflow spectrum is domi¬ 
nated by neutral helium lines and appears to follow preferred 
paths corotating with the binary system. Since the line flux 
along these paths is sporadic, further study is needed to de¬ 
termine the persistent patterns, but emission along a pulsar 
wind-driven spiral outflow provides a viable model. Select¬ 
ing against companion flares and spectrally excluding the out¬ 
flow line flux, we have made improved measurements of the 
companion photosphere’s heated light curve and color varia¬ 
tion. These data strengthen the conclusion that simple radia¬ 
tive heating models are not adequate to describe strongly in¬ 
teracting black-widow binaries like J1311. Instead, we posit 
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Fig . 11.— Confidence regions for two-parameter ELC model fits, holding 
Lfi and Tjv fixed near their maximum values and varying the assumed spot 
deficit from l.Ox (no spot) to 0.5x (deep spot, excluded by the lack of an 
eclipse). The best-fit point and the x^/DoF are marked for each deficit. The 
uncertainty ellipses are for Xmin4(x^T)oF) and Xmin+3^ (x^dDoF). Lines 
of constant Mi (solid, dashed) and M 2 (dotted) are drawn. The inclination 
and mass ratio, and hence the infened pulsar mass, depend dramatically on 
the heating pattern. 

that the heating flux is reprocessed in an intrabinary shock 
whose illumination of the companion surface preferentially 
heats away from FI, toward the Day-Night terminator. If this 
deep heating is mediated by shock-accelerated particles, co¬ 
herent magnetic fields near FI can also help redirect the flux 
and decrease radiation from the FI point. 

The net effect of asymmetric heating is to prefer smaller 
ATcor and larger sin i, but the precise values of these critical 
parameters are very sensitive to the heating pattern. There 
are two important paths to improving their determination. ( 1 ) 
Improved, more physical, modeling can better exploit exist¬ 
ing measurements. However, statistically acceptable will 
probably require excellent exclusion of flare perturbations, 
even if the heating model is physically correct. (2) Obser- 
















10 


Romani et al. 


vationally, one would like to directly measure the heating pat¬ 
tern across the companion Day side. Improved spectroscopy 
remains the key. While well-measured colors, especially near 
minimum {(pB = 0.75), can help constrain the fits, we have 
seen in this study that these can best be obtained by synthe¬ 
sizing colors from spectra, after excluding flares and the inter¬ 
mittent wind emission lines. Of course, with sufficient S/N, 
the spectra themselves will help determine the distribu¬ 
tion as a function of phase, by measuring relative absorption¬ 
line strengths for species with differing excitation. Even more 
directly, higher spectral resolution with good S/N can measure 
rotational broadening of the companion’s absorption lines and 
its variation with orbital phase, thereby providing perhaps the 
cleanest constraints on sin/. 

Until such improvements are made, we cannot exclude any 
mass in the range 1.8-2.7 Mq for PSR J1311-3430. This 
range of values still allows a wide range of system evolu¬ 
tions and may, or may not, constrain the equation of state at 


Abdo, A. A., et al. 2009b, Science, 325, 840 
Abdo, A. A., et al. 2010, ApJS, 187, 460 

Benvenuto, O. G., De Vito, M. A., & Horvath, L. E. 2012, Apt, 753, L33 
Benvenuto, O. G., De Vito, M. A., & Horvath, L. E. 2013, MNRAS, 
submitted. 

Breton, R. P., et al. 2013, ApJ, 769, 108 

Castelli, F. & Kurucz, R. L. 2004, arXiv:astro-ph/0405087 

Covey, K. R., et al. 2007, AJ, 134, 2398 

Demorest, P, et al. 2010, Nature, 467, 1081 

Faber, S. M., et al. 2003, SPIE, 4841, 1657 

Fichtel, K., et al. 1994, ApJS, 94, 551 

Filippenko, A. V. 1982, PASP, 94, 715 

Greiner, J., et al. 2008, PASP, 120, 405 

Hauschildt, P. H., et al. 1999, ApJ, 525, 871 

Hook, I., et al. 2005, PASP, 116, 425 

Jeffery, C. S., Woolf, V. M., & Pollacco, D. L. 2001, A&A, 376, 497 

Kataoka, J., et al. 2012, ApJ, 757, 176 

Kurtz, M. J. & Mink, D. J. 1998, PASP, 110, 934 


high densities. Thus, more work is needed to fully exploit the 
promise of this unusual system. 


We thank Kelsey Clubb, Ori Fox, and Melissa Graham for 
assistance with some of the Keck observations, as well as Ger¬ 
man Gimeno for help with the Gemini campaign. We also 
thank the referee for many detailed comments. This work 
was partially supported by NASA grants NNXl 1A044G and 
NNX12A068G. A.V.F. and S.B.C. were supported by Gary 
and Cynthia Bengier, the Richard and Rhoda Goldman Fund, 
the Christopher R. Redlich Fund, the TABASGO Foundation, 
and NSF grant AST-1211916. Some of the data presented 
herein were obtained at the W. M. Keck Observatory, which is 
operated as a scientific partnership among the California Insti¬ 
tute of Technology, the University of California, and NASA; 
the Observatory was made possible by the generous financial 
support of the W. M. Keck Foundation. 


Lattimer, J. M., & Prakash, M. 2007, PhR, 442, 109L 

Lattimer, J. M., & Schutz, B.F. 2005, ApJ, 629, 797 

Nakamura, Y, & Kitamura, M. 1992, ApSS, 191, 267 

Nolan, P. L., et al. 2012, ApJS, 199, 31 

Oke, J. B., et al. 1995, PASP, 107, 375 

Orosz, J. A., & Hauschildt, P. H. 2000, AA, 364, 265 

Pletsch, H. J., et al. 2012, Science, 338, 1314 

Rasio, F. A., Shapiro, S. L., & Teukolsky, S. A. 1989, ApJ, 342, 934 

Romani, R. W. 2012, ApJ, 754, L25 

Romani, R. W., & Shaw, M. S. 2011, ApJ, 743, L26 

Romani, R. W., et al. 2012, ApJ, 760, L36 (R12) 

Romani, R. W., et al. 2014, ApJ, 793, L20 
Schlally, E. F, & Finkbeiner, D. P. 2011, ApJ, 737, 103 
Schroeder, J., & Halpem, J. 2014, ApJ, 793, 78 
van Haften, L. M., et al. 2012, AA, 537, 104 


